Further explorations of Skyrme-Hartree-Fock-Bogoliubov mass formulas. 

II: Role of the effective mass. 
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We have constructed four new complete mass tables, referred to as HFB-4 to HFB-7, each one 
including all the 9200 nuclei lying between the two drip lines over the range of Z and N > 8 and 
Z < 120. HFB-4 and HFB-5 have the isoscalar effective mass Mg constrained to the value 0.92M, 
with the former having a density-independent pairing, and the latter a density-dependent pairing. 
HFB-6 and HFB-7 are similar, except that Mg is constrained to 0.8M. The rms errors of the mass- 
data fits are 0.680, 0.675, 0.686, and 0.676 MeV, respectively, almost as good as for the HFB-2 mass 
formula, for which M* was unconstrained. However, as usual, the single-particle spectra depend 
significantly on M * . This decoupling of the mass fits from the fits to the single-particle spectra has 
been achieved only by making the cutoff parameter of the 5-function pairing force a free parameter. 
An improved treatment of the center-of-mass correction was adopted, but although this makes a 
difference to individual nuclei it does not reduce the overall rms error of the fit. The extrapolations 
of all four new mass formulas out to the drip lines are essentially the same as for the original HFB-2 
mass formula. 

PACS numbers: 21.10.Dr,21.30.-x,21.60.Jz 
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I. INTRODUCTION 

In the last few years it has become possible to construct 
complete mass tables by the Hartree-Fock (HF) method 
PI 0, LJ LM > with the parameters of the underlying forces 
being fitted to essentially all of the available mass data. 
These HF calculations are based on conventional Skyrme 
forces of the form 



stronger for an odd number of nucleons (V^ q ) than for an 
even number {V^~ q ), i.e., the pairing force between neu- 
trons, for example, depends on whether TV is even or odd. 
The two HFB mass formulas 0, E| add to the energy 
corresponding to the above force (and to the kinetic en- 
ergy and Coulomb energy including the exchange term in 
Slater approximation) a phenomenological Wigner term 
of the form 
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and a 5-function pairing force acting between like nu- 
cleons treated either in the full Bogoliubov framework 
(HFB)jlliEl, or the BCS approximation thereto (HF- 
BCS) a. 
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where p = p(r) is the local density, and po is its equilib- 
rium value in symmetric infinite nuclear matter (INM). 
Only in the most recent paper Q was the possibility of 
a density-dependent pairing force admitted; in the first 
three llly Ll we had r\ = 0. However, in all four HF mass 
formulas the strength parameter V^ q was allowed to be 
different for neutrons and protons, and also to be slightly 
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a somewhat simpler Wigner term was used in the HFBCS 
mass formula pj and the first HFB mass formula Q- 

An important question concerns the cutoff to be ap- 
plied to the (5-function pairing force: both BCS and 
Bogoliubov calculations diverge if the space of single- 
particle (s.p.) states over which such a pairing force is 
allowed to act is not truncated. However, making such a 
cutoff is not simply a computational device but is rather 
a vital part of the physics, pairing being essentially a 
finite-range phenomenon. To represent such an interac- 
tion by a (5-function force is thus legitimate only to the 
extent that all high-lying excitations are suppressed, al- 
though how exactly the truncation of the pairing space 
should be made will depend on the precise nature of the 
real, finite-range pairing force. It was precisely our igno- 
rance on this latter point that allowed us in Ref. [3| to 
exploit the cutoff as a new degree of freedom: we found 
an optimal mass fit with the spectrum of s.p. states £j 
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confined to lie in the range 



Ef - £A < £i < E F + £ A , 



(4) 



where Ep is the Fermi energy of the nucleus in question, 
and £a is a free parameter. We shall adopt the same 
parametrization in the present paper. 

The two most recent of our mass formulas 0, had 
their forces, labeled BSk2 and BSk3, respectively, fitted 
to the 2135 nuclei with Z, N > 8 whose masses have been 
measured and compiled in the 2001 Atomic Mass Eval- 
uation (AME) of Audi and Wapstra 0. The essential 
difference between these two forces is that the pairing is 
density-independent in the case of BSk2 0, while having 
a density dependence of the form @, with the param- 
eters T) and a taking the values given by Garrido at al. 
0, in the case of BSk3 J] (no other choice for rj and a 
leads to a significant improvement). The rms errors of 
these fits are 0.674 MeV and 0.656 MeV 0, respec- 
tively; the slight superiority of the latter is too insignifi- 
cant to imply that the mass data require the pairing to 
be density-dependent, and the most that one can say is 
that a density dependence of the form of Ref. is not 
ruled out. (On the other hand, the simple model of r\ 
= 1, corresponding to vanishing pairing in the nuclear 
interior, is quite incompatible with the mass data.) Us- 
ing these two forces, complete mass tables, referred to as 
HFB-2 and HFB-3 , respectively, were constructed, 
including all the 9200 nuclei lying between the two drip 
lines over the range of Z and N > 8 and Z < 120. 

This paper is the second in a series of studies of pos- 
sible modifications to the original HFB-2 calculation 0, 
with respect to both the force model and the method of 
calculation. Our motivation for making such modifica- 
tions is mainly astrophysical: see Section I of paper I of 
this series , in which we dealt with the question of the 
density dependence of the pairing force. The main pur- 
pose of the present paper is to examine the role of the 
effective nucleon mass in Skyrme-HFB mass formulas. 

We begin by recalling that the HF equation, i.e., the 
equation determining the s.p. states, has for the Skyrme 
forces and our choice of symmetries the particularly 
simple form 
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All quantities appearing here are defined as in, for exam- 
ple, Ref. 0, but the essential point is that all the non- 
locality is confined to the effective mass M*(r), which is 
given in terms of the Skyrme parameters by 
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the subscript q denoting neutron or proton. At any point 
in the nucleus the two effective masses M* (r) and M* (r) 
are now seen to be determined entirely by the local den- 
sities according to 
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where M* and M* are the so-called isoscalar and isovec- 
tor effective masses, respectively, quantities that are de- 
termined by the Skyrme-force parameters according to 
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The single-particle energies €i >q of even nuclei, obtained 
as eigenvalues of the HF Hamiltonian, Eq. (J5J , are often 
identified with the one-nucleon separation energies into 
or from certain low-lying excited states in adjacent odd- A 
nuclei, see 8] and references therein. It is known that to 
have the same density of s.p. levels e^q in the vicinity of 
the Fermi level as observed in experiment for heavy and 
intermediate-mass nuclei, one must have MT/M equal 
to, or close to, 1.0 at saturation density po jjj, 0] (we 
see from Eq. Q that the isovector effective mass M v 
will have little influence on s.p. energies of nuclei that 
are relatively close to the stability line). On the other 
hand, INM calculations with forces that are realistic in 
the sense that they fit the two- and three-nucleon data 
(and therefore require an explicit treatment of the short- 
range correlations that are built in an effective way into 
the forces for HFB calculations) indicate that M* JM 
lies in the range 0.6-0.9 for p = p [H E E Q 111 . 
Rough experimental confirmation that M* is indeed sig- 
nificantly smaller than M came first from measurements 
of the deepest s.p. states in light nuclei (the deepest 
s.p. states of heavier nuclei have not been measured); see 
Refs. 0,0,0 for theoretical discussions. More precise 
empirical information comes from analyses of the giant 
isoscalar quadrupole resonance, which lead to a value of 
around 0.8M for M* at p — p , according to Ref. 0. 

Actually, there is no contradiction between these two 
sets of values of M*/M, since Refs. 0,12^ have shown 
that in finite nuclei one can obtain reasonable s.p. level 
densities near the Fermi level with the INM values of 
M*/M, i.e., of 0.6-0.9, provided one takes into account 
the coupling between s.p. excitation modes and surface- 
vibration RPA modes. Since the good agreement with 
measured s.p. level densities found in Ref. 0] was ob- 
tained without making these corrections it must be sup- 
posed that the resulting error is being compensated by 
the higher value of M*/M, i.e., M*/M ~ 1.0, which 
may thus be regarded as a phenomenological value that 
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leads to good agreement with measured s.p. energies 
in straightforward HF, or other mean-field calculations, 
without any of the complications of Refs. 0, ■ 

Now the fact that in all our previous mass fits 0, 0, 
0, where no constraints were placed on the effective 
mass, we found M*/M ~ 1.0 at the density p = po sug- 
gests that obtaining a correct s.p. spectrum in the vicinity 
of the Fermi level is a necessary condition for an optimum 
mass fit. This conclusion tends to be confirmed by the 
occupation-number representation of the Strutinsky the- 
orem, which approximates the total energy of the nucleus 
as 

E-E + ^eidm , (9) 

i 

where E is a smoothed, average value of E, while 8 rii — 
rii — Hi, in which rii is the actual occupation number of 
the s.p. state i and n, is an average occupation number, 
given by. for example, Eq. (IV. 18) of the review of Brack 
et al. [23j. But in Section IV. 6 of this same paper |23( it is 
shown that Srii is non-vanishing only for s.p. states lying 
within about 20 A -1 / 3 MeV of the Fermi level. We may 
expect that it will be difficult to obtain correct masses 
if the s.p. spectrum over this interval is not reproduced, 
and it is quite comprehensible that the optimal mass fits 
published so far require that M*/M take a value close to 
1.0 at the density p — po. Indeed, with a Skyrme force 
of the form Q and density-independent (5-function pair- 
ing forces treated within the BCS approximation, Farine 
et al. 24] found that if M*/M was constrained to be 
equal to 0.8 then in fitting 416 quasi-spherical nuclei it 
was impossible to reduce the rms error below 1.141 MeV 
(parameter set MSk5*). On the other hand, when the 
constraint on M*/M was released the rms error for the 
same data set fell to 0.709 MeV, M*/M rising to 1.05 
(parameter set MSk5). Moreover, the Skyrme force SLy4 
|25j with M*/M = 0.7 does not reproduce very well the 
masses of the few open-shell nuclei considered (see Figs. 
l-4ofRef. |2U). 

Nevertheless, there remains some room for maneuver in 
Eq. 0, since shifts in the s.p. energies ej resulting from 
a change in the effective mass could in principle be com- 
pensated by appropriate changes in the 5 rii, provided full 
use was made of all the degrees of freedom in the force, 
with particular emphasis on those that have not hitherto 
been exploited; of especial interest in this respect is the 
pairing cutoff, which, for example, was fixed at the same 
value in the MSk5 and MSk5* fits (41 A" 1 / 3 MeV into 
the continuum). In the present paper we pursue these 
possibilities with the aim of seeing to what extent it is 
possible to maintain quality mass fits with lower values 
of M*/M, more appropriate to INM. Part of our own 
interest in this question lies in our search for a unique 
effective interaction suitable for the determination of an 
equation of state describing the formation of nuclear mat- 
ter from isolated finite nuclei that occurs during stellar 
collapse 01131 ■ A lower value of M*/M, corresponding 
to INM, will certainly be appropriate in the final stage; 



the question here is to see to what extent such a choice 
of M*/M can suitably describe the isolated nuclei pre- 
vailing at the beginning of the collapse. (Similar consid- 
erations will arise in the more or less inverse sequence 
of events traced out during the neutron-matter decom- 
pression that occurs, for example, in the aftermath of 
neutron-star mergers.) 

In this paper we also replace our original approximate 
correction for the centre-of-mass motion by a much im- 
proved treatment, and we describe this first, in Section II, 
leaving until Section III the discussion of the role of the 
effective mass. There, two different values of M* at the 
density p — po, 0.92M and 0.8M, will be considered, and 
four new mass tables, HFB-4, HFB-5, HFB-6, and HFB- 
7, generated, each value of M* being calculated with and 
without a density dependence in the pairing force. 



II. CENTER-OF-MASS CORRECTION 

Mean-field approaches like HF or HFB, which establish 
an intrinsic frame of the nucleus, break several symme- 
tries of the Hamiltonian and the wave function in the 
laboratory frame 0,|22j ■ For example, finite nuclei break 
translational invariance, deformed nuclei are not rota- 
tional invariant, and the HFB approach breaks particle- 
number symmetry. Doing so adds desired correlations to 
the modeling - as multi-particle-multi-hole states associ- 
ated with deformation, pairing, etc. - but at the same 
time gives rise to an admixture of excited states to the 
calculated ground state. Their spurious contribution to 
the total mass changes with nucleon numbers and defor- 
mation. The energy of those spurious modes, which are 
not explicitly removed from the calculated ground state, 
will be simulated by the Skyrme force through the param- 
eter fit, which might spoil the properties of the resulting 
forces [29|. 

A rigorous way to restore the broken symmetries is pro- 
jection on exact quantum numbers, but this would be too 
time-consuming to be used for the large-scale mass fits 
performed here. A simpler procedure is to estimate the 
contribution to the binding energy in a suitable approxi- 
mation, and to add the resulting corrections to the calcu- 
lated masses; such a procedure has already been adopted 
by many workers for the center-of-mass (cm) correction, 
the rotational correction, and the Lipkin-Nogami correc- 
tion to the pairing energy. In our own HF mass formulas 
IH IE H we have used such a procedure for the cm. 
and rotational corrections, as described in Ref. |7j . In the 
present paper we improve our treatment of the cm cor- 
rection, as discussed below in this section, but otherwise 
the Skyrmc-HFB formalism used here is essentially as de- 
scribed in detail in our first HFB paper 2]. In particular, 
we have not yet made any correction for particle-number 
fluctuations; this will be the topic of a forthcoming paper. 

The HFB ground state is not an eigenstate of the 
total momentum operator. Thus, although the expec- 
tation value of the momentum operator P = pi 
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FIG. 1: Comparison for the 2135 nuclei included in the 2001 
mass compilation of |5j of the binding energies obtained with 
the improved centre-of-mass correction E cm with those ob- 
tained with the approximation of [2gj| E app . Both calculations 
use the BSk2 Skyrme force and are made in the spherical ap- 
proximation. 
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FIG. 2: Differences between the HFB-4' and HFB-4 masses 
as a function of the neutron number N for all nuclei with 
Z, N > 8 lying between the proton and neutron driplines up 
to Z = 120. 



in the cm frame (HFB|P|HFB) vanishes, its dispersion 
(HFB|P 2 |HFB) does not. Gaussian overlap approxima- 
tion to exact momentum projection gives for the spurious 
cm energy 

Ecm = 2Abi (HFB| ^ 2|HFB> ' (10) 

which has to be subtracted from the calculated total en- 
ergy. To avoid the time-consuming evaluation of the non- 
diagonal terms jij ■ pj, i ^ j, in the past we have always 
adopted the approximation of Butler et al. psj . which 
takes explicit account only of the diagonal terms, thus 

£ cm ^^(HFB|^p 2 |HFB) , (11) 

i 

where f(A) is a simple function that makes this expres- 
sion exact in the case of pure oscillator s.p. states. From 
now on, we evaluate the cm correction according to Eq. 
I) l()[l. doing so, however, perturbatively. That is, both the 
diagonal and off-diagonal terms of Eq. H10|) are included 
only in the calculation of the converged total energy, not 
in the variational equation that leads to the mean field 
in the HF equation JSJ. 

The effect of this improved treatment is seen in Fig. ^ 
where for each of the 2135 nuclei of known mass we show 
the difference between the total energy calculated with 
our new method, and that calculated with the approxi- 
mation of Ref. 28]. The force used for this comparison 
is BSk2, and for simplicity we assume a spherical config- 
uration for all nuclei. The differences are largest for light 
nuclei, for which they can reach 1.5 MeV; strong shell 
effects will be noticed. 

Of course, in fitting force BSk2 to the mass data these 
errors in the approximation of Ref. psj were absorbed 
to some extent into the force parameters. It would be 
interesting to refit the BSk2 force to the same data us- 
ing the improved cm correction, and compare the new 



with the original BSk2 fit. For practical reasons, wc 
make this comparison on the basis of the force BSk4, 
described in the next section. This force (like the forces 
BSk5, BSk6 and BSk7) is calculated with the improved 
cm correction |29j , so we refit the force BSk4 to the data 
using the approximate cm correction of Ref. |28|. defin- 
ing thereby the force BSk4'. We find that the change 
in the rms error of the fit is negligible, going from 0.680 
to 0.681 MeV. The mass differences between the corre- 
sponding BSk4' and BSk4 predictions are compared in 
Fig. |21 for all nuclei with Z, N > 8 lying between the pro- 
ton and neutron driplines up to Z — 120. In contrast to 
Fig. ^ deformation effects are taken into account consis- 
tently. The difference in the shell structure between the 
improved treatment of the cm correction pflj and the ap- 
proximation of Ref. j2^, as seen in Fig.^ is still present 
after renormalizing the Skyrme forces on experimental 
masses. The mass differences remain however smaller 
than 2 MeV, even for exotic neutron-rich or superheavy 
nuclei. 



III. CHOICE OF THE ISOSCALAR EFFECTIVE 
MASS, AND THE NEW HFB MASS TABLES 

Most INM calculations of M* at the density p = po 
give a value of around 0.8M, the most recent such calcu- 
lation being that of Ref. [15| . On the other hand, using 
the so-called extended Brueckner-Hartree-Fock method 
with realistic nucleonic forces, Ref. |b| finds 0.92A/. 
Rather than attempt to decide between these two values 
we shall here consider both of them, with the value of 
0.92Ajf being imposed on parameter sets BSk4 and BSk5 
(mass formulas HFB-4 and HFB-5, respectively), and the 
value of 0.8M being imposed on parameter sets BSk6 and 
BSk7 (mass formulas HFB-6 and HFB-7, respectively). 
The pairing force in both BSk4 and BSk6 is supposed 
to be density-independent, while that of BSk5 and BSk7 
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TABLE I: Rms (<r) and mean (e) errors (in MeV) in the predictions of masses M obtained with the BSk2-7 forces. Also given 
are the model standard deviation and the model mean error (see text for more details) on the set of all 2135 measured masses, 
on the 928 stable and neutron-rich nuclei, and on the 1207 proton-rich nuclei. The last two lines correspond to the rms and 
mean errors (in fm) for the predictions of the 523 measured charge radii (r c ). 





BSk2 


BSk3 


BSk4 


BSk5 


BSk6 


BSk7 


cr(M) (2135 nuclei) 


0.674 


0.656 


0.680 


0.675 


0.686 


0.676 


e(M) (2135 nuclei) 


0.000 


-0.006 


-0.115 


-0.005 


-0.013 


-0.004 


cr mo d(M) (2135 nuclei) 


0.660 


0.639 


0.661 


0.655 


0.666 


0.658 


£mod(M) (2135 nuclei) 


-0.007 


-0.015 


0.106 


-0.006 


0.013 


0.026 


(T mo d(M) (928 nuclei) 


0.709 


0.709 


0.678 


0.685 


0.713 


0.707 


e mod (M)(928 nuclei) 


-0.168 


-0.118 


0.001 


-0.122 


0.093 


-0.085 


o mod {M) (1207 nuclei) 


0.620 


0.581 


0.649 


0.631 


0.629 


0.618 


e mod (M)(12Q7 nuclei) 


0.113 


0.062 


0.186 


0.080 


0.092 


0.109 


a(r c ) (523 nuclei) 


0.0282 


0.0291 


0.0282 


0.0270 


0.0262 


0.0260 


e(r c ) (523 nuclei) 


0.0138 


0.0161 


-0.0115 


0.0104 


0.0028 


-0.0030 



TABLE II: Skyrme-force and pairing-force parameters of BSk2-BSk7 





BSk2 


BSk3 


BSk4 


BSk5 


BSk6 


BSk7 


t [McV fm 3 ] 


-1790.6248 


-1755.1297 


-1776.9376 


-1778.8934 


-2043.3174 


-2044.2484 


t\ [McV fm 5 ] 


260.996 


233.262 


306.884 


312.727 


382.127 


385.973 


t 2 [MeV fm 5 ] 


-147.167 


-135.284 


-105.670 


-102.883 


-173.879 


-131.525 


t 3 [McV fm 3+3T ] 


13215.1 


13543.2 


12302.1 


12318.37 


12511.7 


12518.8 


X 


0.498986 


0.476585 


0.542594 


0.444510 


0.735859 


0.729193 


Xl 


-0.089752 


-0.032567 


-0.535165 


-0.488716 


-0.799153 


-0.932335 


X 2 


0.224411 


0.470393 


0.494738 


0.584590 


-0.358983 


-0.050127 


X3 


0.515675 


0.422501 


0.759028 


0.569304 


1.234779 


.236280 


Wo [McV fm 5 ] 


119.047 


116.07 


129.50 


130.70 


142.38 


146.93 


7 


0.343295 


0.361194 


1/3 


1/3 


1/4 


1/4 


V„ [McV fm 3 ] 


-238 


-359 


-273 


-429 


-321 


-505 


V~ [McV fm 3 ] 


-265 


-407 


-289 


-463 


-325 


-514 


V+ [McV fm 3 ] 


-247 


-365 


-285 


-447 


-338 


-531 


V~ [McV fm 3 ] 


-278 


-413 


-302 


-483 


-341 


-541 


V 





0.45 





0.45 





0.45 


a 





0.47 





0.47 





0.47 


£A[McV] 


15 


14 


16 


16 


17 


17 


Vw [MeV] 


-2.05 


-2.05 


-1.72 


-1.96 


1.76 


1.86 


A 


485 


460 


740 


480 


700 


720 


V w [McV] 


0.70 


0.54 


0.54 


0.50 


0.58 


0.54 


A 


28 


30 


30 


30 


28 


28 



is taken to have a density dependence of the form J3J, 
with the parameters r\ and a taking the values given by 
Garrido at al. [||, as with the BSk3 force In all four 
cases the isovector mass at p — po is left unconstrained 
in the fits. 

We fit all four of these forces to the same data set 
as the one to which BSk2 3] and BSk3 were fitted, 
i.e., the 2135 nuclei with Z, N > 8 whose masses have 
been measured and compiled in the 2001 AME As 
in the fits of BSk2 and BSk3, we impose a lower limit 
on the INM symmetry coefficient J, in order to prevent 
the collapse of neutron matter at nuclear densities, as 



required by the observed stability of neutron stars. In 
the case of BSk5, we use the INM calculation of 0] 
to constrain not only the effective mass to 0.92M, but 
also the symmetry coefficient to J = 28.7 MeV. Also, we 
set the equilibrium density of symmetric INM at po = 
0.1575 fm -3 , it having been found that this ensures very 
good predictions for rms charge radii along with near- 
optimal mass fits. 

From the first line of Table [Q we see that while forces 
BSk2 and BSk3, in which M* is unconstrained, still give 
the best fits, all four of the new forces, BSk4-7, give fits 
that are almost as good, the deterioration on constraining 
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FIG. 3: Differences between experimental and calculated 
mass excesses as a function of the neutron number N for the 
HFB-4 to HFB-7 mass tables. 



M* being quite negligible. Fig. [^displays the deviations 
of the calculated masses from the experimental values. 

Table [I] also has entries for a "model standard devia- 
tion" <7 mo d and a "model mean error" l mo d- These quan- 
tities have been introduced [3(], |^| as an improved mea- 
sure of the validity of the physical model that is being 
fitted to the data, the standard rms error a suffering from 
the defect that it does not take account of the experimen- 
tal errors of the individual mass measurements, as given 
by Audi and Wapstra @, |3^] . The standard rms error 
is a legitimate measure when the experimental errors are 
small compared to the rms error itself, but some of the 
most recent measurements of nuclei far from the stabil- 
ity line have errors in excess of I MeV @ . In such cases, 
o'mod and lmod give a better assessment of the validity of 
a given mass formula, since they weight each data point 
in terms of its experimental error, following a procedure 
based on the method of maximum likelihood [M |3l|. 
The definition of a mo d that we adopt is that of Eqns. 



(42) and (43) of Ref. [3(j (which writes cr mo rf as (T^ h ), 
while our definition of e mo ,j is that of Eqns. (9) and (10) 
of Ref. lU (which writes e mo d as fJ^ h ); for a discussion of 
the relationship between the ways in which the two pa- 
pers y(J and [3l| treat model errors see Appendix B of 



Ref. 



I . We show Omod and e mo w for the full data set of 



2135 masses to which the fit was made, as well as to the 
two subsets of 1207 proton-rich nuclei and 928 stable and 
neutron-rich nuclei. None of these model errors suggests 
that any particular force is significantly better or worse 
than any of the others, as far as masses are concerned. 
The fact that mass fits with values of M* /M constrained 
to be equal to 0.8 can be obtained which are almost as 
good as those for which M*/M is unconstrained (always 
emerging with a value close to 1.0) is, we have found, es- 
sentially a consequence of our exploitation of the degree 
of freedom associated with the pairing cutoff, a reduc- 
tion in M* jM being almost completely compensated by 
an increase in the parameter ea appearing in Eq. |@J (see 
Table CU). 

We also show in Tabled the rms and mean deviations 
between our calculated and experimental charge radii for 
the 523 nuclei listed in the 1994 compilation |34[ (for 
more details on the HFB derivation of the charge radii, 
see Ref. [35| ) . The overall agreement with experiment is 
seen to be excellent. However, none of the forces is able 
to completely reproduce the much discussed kink in the 
Pb isotope chain at N = 126. 

The parameters of the forces BSk2-7 are given in Ta- 
ble [HJ while the corresponding macroscopic parameters, 
i.e., the parameters relating to INM and SINM (semi- 
infinite nuclear matter) calculated for all these forces, 
are shown in Table ITTT1 The quantities appearing in this 
latter table that have not yet been defined are: a v , the 
energy per nucleon at equilibrium in symmetric INM; 
K v , the INM incompressibility; Go and G' , the Lan- 
dau parameters defined in Ref. |3(|; pfrmgi the density 
at which neutron matter flips over into a ferromagnetic 
state that has no energy minimum and would collapse 
indefinitely |37l ]: a s f, the surface coefficient; and Q, the 
surface-stiffness coefficient |3g. It will be recalled that 



TABLE III: Macroscopic parameters of the forces BSk2-BSk7 





BSk2 


BSk3 


BSk4 


BSk5 


BSk6 


BSk7 


a v [McV] 


-15.794 


-15.804 


-15.773 


-15.802 


-15.749 


-15.760 


po [fm~ 3 ] 


0.1575 


0.1575 


0.1575 


0.1575 


0.1575 


0.1575 


J [McV] 


28.0 


27.9 


28.0 


28.7 


28.0 


28.0 


m;/m 


1.04 


1.12 


0.92 


0.92 


0.80 


0.80 


M*/M 


0.86 


0.89 


0.85 


0.84 


0.86 


0.87 


K v [MeV] 


233.6 


234.8 


236.8 


237.2 


229.1 


229.3 


Go 


-0.705 


-0.994 


-0.478 


-0.579 


0.065 


-0.101 


Go 


0.446 


0.497 


0.457 


0.454 


0.312 


0.356 


Pfrmg/po 


1.1 


1.00 


1.30 


1.20 


1.81 


1.62 


a s f [McV] 


17.5 


17.5 


17.3 


17.5 


17.2 


17.3 


Q [McV] 


68 


67 


76 


52 


83 


80 
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in all cases the values of po and J were imposed, as were 
the values of M*/M. On the other hand, M*/M is un- 
constrained, but it is gratifying to note that for all fits 
its value comes out to be quite close to the value of 0.83 
that we infer from the INM calculations of Zuo et al. 0] 
(see especially their Fig. 9). Nevertheless, a word of cau- 
tion is necessary here, since it is known from Ref. [39| 
that the rms error of the mass fit varies only slowly with 
M*. As for K v , all our forces give values falling within 
the experimental rang e of 225-240 MeV established by 
Youngblood et al. |4(| • All our forces likewise satisfy the 
condition Go and Gn > — 1 for stability against spin and 
spin-isospin flips |4l| at saturation density. Finally, it 
will be seen that reducing M* seems to assure greater 
stability to neutron matter against ferromagnetic flips. 
For the values of pf rm g, Gq and G' given in Table ITTT1 it 
is assumed that the effective spin-spin interaction is ob- 
tained from the exchange terms of the two-body Skyrme 
force, eq. . There is another possible view of the effec- 
tive Skyrme interaction, see H, ar± d references given 
therein, which leaves the coupling constants that deter- 
mine Pfrmg, Gq and G' as additional free parameters not 
constrained by our current mass fit. 

Neutron Matter. All our present forces have a value of 
J ~ 28 MeV (except BSk5 with J = 28.7 MeV) which 
was actually set as a lower limit in the search on the 
Skyrme parameters. It is possible that a slightly better 
mass fit could have been obtained with a value somewhat 
closer to 27.5 MeV, but this might have engendered an 
unphysical collapse of neutron matter at nuclear densi- 
ties. The situation in neutron matter for our forces is 
as shown in Fig. the solid curve labeled FP shows 
the results of Friedman and Pandharipande 0] f° r the 
realistic force Argonne V14 + TNI, containing two- and 
three- nucleon terms. More recent realistic calculations of 
neutron matter 0> 19 S^ ve similar results up to nu- 
clear densities. At low densities, all forces give the same 
predictions. At densities higher than the nuclear den- 
sity, where the validity of the non-relativistic approach 
followed here could admittedly be questioned, the new 




FIG. 4: Energy-density curves of neutron matter for the forces 
of this paper, and for the calculations of Ref. 1121 (FP). 



TABLE IV: Single-particle proton levels in 20S Pb (MeV). Ex- 
perimental values are taken from Ref. [r?|. The asterisk de- 
notes the Fermi level. The quantity A p is the interval between 
the centroids of the lg and 2/ doublets. 



Level 


BSk2 BSk3 BSk6 MSk5 MSk5* SLy4 


Expt. 


lSl/2 


-33.1 


-31.4 


-38.9 


-31.8 


-40.5 


-44.0 


- 


lff9/2 


-14.8 


-14.4 


-16.2 


-14.6 


-16.7 


-17.7 


-15.4 


107/2 


-11.4 


-11.1 


-12.6 


-11.3 


-13.1 


-13.5 


-11.4 


2d 5/2 


-9.8 


-9.6 


-10.4 


-9.7 


-10.5 


-11.5 


-9.7 


lft-11/2 


-8.8 


-8.7 


-9.0 


-8.7 


-9.2 


-9.7 


-9.4 


2^3/2 


-8.2 


-8.1 


-8.5 


-8.2 


-8.8 


-9.6 


-8.4 


3S1/2* 


-7.6 


-7.5 


-7.9 


-7.6 


-7.9 


-8.8 


-8.0 


1^9/2 


-3.9 


-3.9 


-3.8 


-4.0 


-4.1 


-3.8 


-3.8 


2/7/2 


-3.1 


-3.3 


-2.6 


-3.3 


-2.3 


-2.9 


-2.9 


1*13/2 


-2.2 


-2.4 


-1.5 


-2.4 


-1.3 


-1.5 


-2.2 


3^3/2 


-0.3 


-0.6 


0.6 


-0.6 


0.9 


0.4 


-1.0 


2/5/2 


-0.9 


-1.1 


0.0 


-1.1 


0.0 


-0.4 


-0.5 


A p 


11.0 


10.5 


13.1 


10.7 


13.9 


14.0 


11.7 



TABLE V: Single-particle neutron levels in 208 Pb (MeV). Ex- 
perimental values are taken from Ref. [TH ]. The asterisk de- 
notes the Fermi level. The quantity A n is the interval between 
the centroids of the 2/ and 3d doublets. 



Level 


BSk2 BSk3 BSk6 MSk5 MSk5* SLy4 


Expt. 


lsi/2 


-39.5 


-36.9 


-51.5 


-40.7 


-49.9 


-57.9 




1^9/2 


-10.6 


-10.2 


-12.4 


-10.8 


-12.3 


-12.5 


-10.9 


2/7/2 


-10.8 


-10.6 


-11.7 


-10.8 


-11.5 


-12.0 


-9.7 


1*13/2 


-9.2 


-9.1 


-9.6 


-9.1 


-9.4 


-9.6 


-9.0 


3P3/2 


-8.4 


-8.3 


-8.9 


-8.5 


-8.7 


-9.2 


-8.3 


2/5/2 


-8.2 


-8.1 


-8.8 


-8.3 


-8.7 


-9.1 


-8.0 


3P1/2* 


-7.5 


-7.4 


-7.8 


-7.6 


-7.7 


-8.1 


-7.4 


2#9/2 


-4.2 


-4.3 


-3.6 


-4.1 


-3.4 


-3.2 


-3.9 


1*11/2 


-2.6 


-2.7 


-2.4 


-2.7 


-2.6 


-1.7 


-3.2 


-1.715/2 


-2.4 


-2.6 


-1.3 


-2.1 


-1.2 


-0.6 


-2.5 


3^5/2 


-1.9 


-2.1 


-1.2 


-1.8 


-1.0 


-0.7 


-2.4 


4si/ 2 


-1.2 


-1.4 


-0.6 


-1.0 


-0.3 


0.0 


-1.9 


2^7/2 


-1.2 


-1.4 


-0.2 


-1.1 


-0.3 


0.0 


-1.5 


3^3/2 


-0.8 


-1.0 


0.0 


-0.7 


0.1 


0.3 


-1.4 


An 


8.2 


7.8 


9.8 


8.3 


9.7 


10.5 


7.2 



forces are found to lead to harder neutron-matter curves 
and avoid the collapse obtained for BSk2. In particular, 
BSk5 with J = 28.7 MeV gives very similar energies per 
nucleon as BSk6 and BSk7 with J = 28 MeV due to 
the much lower value of the aq parameter reached by the 
mass fit (see Table |TT| . 
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Single-particle spectra. Clearly, it is of interest to see 
what happens to the s.p. energies when the effective 
mass is reduced while maintaining a very good fit to 
the masses. In Tables IIVI - IIXI and Figure we show 
the s.p. spectra of 208 Pb, 132 Sn, and 16 for forces BSk2, 
BSk3, and BSk6, along with those of the old forces MSk5, 
MSk5* and SLy4. Comparing BSk2 and BSk3 shows 
that making the pairing density-dependent has little ef- 
fect on the s.p. energies. Also, the spectra of BSk2 and 
MSk5 are very similar, but different from those of BSk6, 
MSk5*, and SLy4, which, however, resemble each other 
quite closely. That is, the effective mass still determines 
the s.p. spectra. But while BSk2 and BSk6 have different 
s.p. spectra they give very similar fits to the mass data. 



TABLE VI: Single-particle proton levels in 132 Sn (MeV). Ex- 
perimental values are taken from Refs. [it! |47|. The asterisk 
denotes the Fermi level. The quantity A p is the interval be- 
tween the centroids of the lg and 2d doublets. 



Level 


BSk2 BSk3 BSk6 MSk5 MSk5* SLy4 


Expt. 




-38.7 


-37.1 


-44.2 


-37.5 


-48.1 


-49.0 




2fl/2 


-16.0 


-15.9 


-16.4 


-15.9 


-18.7 


-17.6 


-16.1 


lff9/2* 


-14.7 


-14.6 


-14.9 


-14.7 


-16.6 


-15.6 


-15.8 


1^7/2 


-9.5 


-9.6 


-9.4 


-9.6 


-11.1 


-9.3 


-9.7 


24/2 


-9.5 


-9.7 


-8.9 


-9.7 


-10.9 


-9.2 


-8.7 


24/2 


-7.6 


-7.9 


-6.7 


-7.9 


-8.5 


-6.9 


-7.2 


3S1/2 


-7.3 


-7.6 


-6.3 


-7.6 


-8.7 


-6.4 




lft-11/2 


-7.1 


-7.3 


-6.2 


-7.4 


-7.5 


-6.2 


-6.8 


A p 


3.7 


3.4 


4.5 


3.7 


4.3 


4.5 


5.1 



On the other hand, the s.p. spectra of BSk6 are very sim- 
ilar to those of MSk5*, but the latter gives a much worse 
mass fit. The only way in which we can reconcile this 
behavior with the Strutinsky theorem in the form is 
to invoke the 8ni quantities: shifts in the Srii compen- 
sate the differences between the s.p. spectra of BSk2 and 
BSk6, and at the same time account for the different mass 
predictions of BSk6 and MSk5*, despite their similar s.p. 
spectra. This interpretation in terms of the occupation 
numbers is strengthened by the fact that the decoupling 



TABLE VII: Single-particle neutron levels in 132 Sn (MeV). 
Experimental values are taken from Ref. [47| . The asterisk 
denotes the Fermi level. The quantity A„ is the interval be- 
tween the centroids of the 2d and 3p doublets. 



Level 


BSk2 BSk3 BSk6 MSk5 MSk5* 


SLy4 


Expt 


lSl/2 


-37.8 


-35.3 


-49.8 


-39.0 


-48.1 


-55.8 




lff7/2 


-9.2 


-8.8 


-11.2 


-9.5 


-11.1 


-11.4 


-9.8 


2d 5/2 


-10.2 


-9.9 


-11.2 


-10.3 


-10.9 


-11.7 


-9.0 


3Si/2 


-8.3 


-8.1 


-8.9 


-8.4 


-8.7 


-9.4 


-7.7 


lft-11/2 


-7.3 


-7.2 


-7.5 


-7.1 


-7.5 


-7.7 


-7.6 


2d 3/2 * 


-7.9 


-7.8 


-8.6 


-8.0 


-8.5 


-9.1 


-7.4 


2/7/2 


-2.9 


-3.1 


-2.3 


-2.8 


-2.2 


-2.0 


-2.4 


3f>3/2 


-1.3 


-1.5 


-0.7 


-1.2 


-0.5 


-0.3 


-1.6 


1^9/2 


0.0 


-0.1 


0.2 


-0.1 


0.1 


0.8 


-0.9 


3P1/2 


-0.7 


-0.9 


0.0 


-0.5 


0.1 


0.3 


-0.8 


2/5/2 


-0.4 


-0.6 


0.4 


-0.3 


0.4 


0.6 


-0.4 


A n 


8.2 


7.8 


9.7 


8.4 


9.6 


10.6 


7.1 
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TABLE VIII: Single-particle proton levels in 16 O (MeV). Ex- 
perimental values are taken from Ref. The asterisk de- 
notes the Fermi level. The quantity A p is the interval between 
the centroid of the lp doublet and the 2s 1 / 2 state. (Errors in 
Table 6a of Ref. :24] for MSk5* corrected.) 



Level 


BSk2 BSk3 BSk6 MSk5 MSk5* SLy4 


Expt. 


lSl/2 


-25.5 -24.5 -29.4 -25.4 


-30.3 


-32.9 


-40 ±8 


lf3/2 


-15.1 -14.8 -15.8 -15.1 


-16.6 


-17.2 


-18.4 


lPl/2* 


-9.5 -9.3 -9.4 -9.6 


-10.6 


-11.1 


-12.1 


lds/2 


-4.3 -4.4 -3.1 -4.3 


-3.6 


-3.4 


-0.6 


2Sl/2 


-1.7 -2.0 -0.4 -1.7 


-0.7 


-1.1 


-0.1 


A p 


11.5 11.0 13.2 11.6 


13.9 


14.1 


16.2 



TABLE IX: Single-particle neutron levels in 16 O (MeV). Ex- 
perimental values are taken from Ref. [Tt1| . The asterisk de- 
notes the Fermi level. The quantity A n is the interval between 
the centroid of the lp doublet and the 2si/2 state. 



Level 


BSk2 BSk3 BSk6 MSk5 MSk5* SLy4 


Expt. 


lSl/2 
lf3/2 
lPl/2* 
1^5/2 
2S1/2 


-28.9 -27.9 -33.0 -28.9 -33.9 -36.7 
-18.5 -18.2 -19.2 -18.5 -20.1 -20.7 
-12.7 -12.6 -12.7 -12.9 -13.9 -14.5 
-7.5 -7.7 -6.3 -7.5 -6.9 -6.6 
-4.8 -5.1 -3.3 -4.8 -3.7 -4.0 


-21.8 
-15.7 
-4.1 
-3.3 


A p 


11.8 11.2 13.7 11.8 14.8 14.6 


16.5 



of the mass fits from the fits to the s.p. spectra is made 
possible only by adjustment of the pairing cutoff. 

As for the agreement with the experimental s.p. spec- 
tra, we see that a value of M* /M close to 1.0 is favored by 
the 208 Pb data, while 16 favors a value of 0.8. The 132 Sn 
data are ambiguous, the neutron spectrum indicating the 
higher value of M* /M, and the proton spectrum the 
lower value. Presumably, if we took into account the cou- 
pling of s.p. excitations and surface- vibration RPA modes 
for the forces with M*/M = 0.8, as in Refs. [UHl], the 
calculated s.p. spectra of heavy nuclei would be in better 
agreement with experiment, but we would then have to 
refit to the mass data. 



Extrapolation to drip lines. With each of the forces 
BSk4-7 determined as described we constructed complete 
mass tables, labeled HFB-4 to HFB-7, respectively, for 
the same nuclei as were included in the HFB-2 and HFB- 
3 tables, i.e., all the 9200 nuclei lying between the two 
drip lines over the range of Z and N > 8 and Z < 120. 
The differences between the HFB-2 masses and the HFB- 
4, 5, 6, and 7 masses are displayed in Fig. [S]as a function 




M(HFB-l) -MfHFB-2) 




50 LOO L50 2O0 250 



FIG. 6: Differences between the HFB-2 and HFB-4 to HFB-7 
masses as a function of the neutron number N for all nuclei 
with Z, N > 8 lying between the proton and neutron driplines 
up to Z = 120. 



of the neutron number N and in Fig. \7\&s a function of 
the neutron separation energy S n , where it will be seen 
that these differences never exceed 6 MeV, and even then 
only as the neutron drip line is approached, no matter 
what the value of Z (different published mass formulas 
all giving very good data fits can differ by up to 15 or 
20 MeV at the drip lines [13 )■ Looking at the first and 
third panels of Fig. shows that as M*/M is reduced 
there is a definite tendency for open-shell nuclei to be 
bound a little more strongly, a trend that becomes more 
conspicuous for the heaviest nuclei. The second and last 
panels of this figure confirm the feature already noted 
in paper I 4J for density-dependent pairing: a similar 
tendency for open-shell nuclei to be more strongly bound, 
especially for the heaviest nuclei. 

Actually, while differences of up to 6 MeV between the 
different mass predictions for nuclei far from the stability 
line may appear to be rather large, of far greater interest 
for practical applications such as to the r-process of nucle- 
osynthesis are differential quantities such as the neutron 
separation energy S n . For these quantities the differences 
between the different predictions are much smaller, as is 
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seen in Figs. |S1 and El where we plot as a function of Z, 
for each of the magic numbers No = 50, 82, 126, and 184, 
respectively, the neutron-shell gaps, defined by 

A(N ) ee S 2 n(Z,N )-S 2 „(Z,N + 2) 

= M(Z, N -2) + M(Z, N + 2)- 2M(Z, 7V Q12) 

calculated with BSk2, BSk4 and BSk6 (S^n denotes the 
2-neutron separation energy). From Figs. |S| and |5| it can 
be seen that the gaps do not depend significantly on the 
effective nucleon mass. The impact of such differences 
on the predicted r-process abundance distributions will 
be studied in more detail in a forthcoming paper. 

For Nq = 50, 82, and 126, Figs. EElshow a strong dis- 
agreement with experiment for the new mass formulas in 
the vicinity of the (semi-) magic proton numbers Z = 40, 
50, and 82. This is related to the problem of "mutually 
supporting magicities" that we have already discussed in 
connection with the HFB-2 mass formula 0. Clearly, 
it has not been solved, either by reducing the effective 
mass, or by introducing the shell-dependent cm correc- 
tion. (We showed in paper I Q that making the pairing 
density-dependent cannot help in this respect, either.) 



IV. CONCLUSIONS 

Fitting Skyrme-type forces to the available mass data 
without any constraint on the effective mass always leads 
to an isoscalar effective mass M* close to the real nucleon 
mass M. However, we have shown that we can reduce 
M*/M to 0.8 without any significant reduction in the 
quality of the mass-data fit, although important changes 
in the s.p. spectra are thereby induced. This decoupling 
of the fit to the mass data from the fit to the s.p. data 
was made possible only by exploiting the pairing-force 
cutoff. 

On this basis we constructed four new complete mass 
tables, referred to as HFB-4 to HFB-7, each one including 
all the 9200 nuclei lying between the two drip lines over 
the range of Z and N > 8 and Z < 120. HFB-4 and 
HFB-5 have M*/M constrained to the value 0.92, with 
the former having a density-independent pairing, and the 
latter a density-dependent pairing. HFB-6 and HFB- 
7 are similar, except that M*/M is constrained to 0.8. 
The mass-data fits are almost as good as those given by 
mass formulas HFB-2 and HFB-3, in which M*/M was 
unconstrained. Actually, in these four new mass formulas 
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FIG. 9: N = 126 (upper panel) and N = 184 (lower panel) 
shell gap as a function of Z. 



we have used an improved treatment of the center-of- 
mass correction [2j| , but although this makes a difference 



to individual nuclei we have shown that the overall rms 
errors would have been essentially the same if we had 
used the same correction as in HFB-2 and HFB-3. 

The extrapolations out to the neutron-drip line of all 
these different mass formulas are essentially equivalent. 
We thus see that the mass predictions required for the 
elucidation of the r-process are beginning to acquire a 
certain stability against changes in the underlying model. 
Nevertheless, it must be remembered that the acquisition 
of new mass data in the regions far from stabilty may well 
necessitate drastic changes to the underlying model. 

Although the forces presented in this paper are equiv- 
alent from the standpoint of nuclear masses, there may 
still be significant differences as far as other quantities 
of astrophysical significance are concerned, e.g., fission 
properties, nuclear level densities, giant isovector dipole 
resonance (GDR), and beta-strength functions. Investi- 
gations along these lines has already begun, and it has 
been shown [43 that the measured positions of the GDR 
strongly favor the Skyrme forces BSk6 and BSk7 with 
their low effective mass of M * — 0.8M ; these calculations 
were made within the HFB plus Quasi-particle Random 
Phase Approximation (QRPA) framework (the second- 
RPA method being applied to estimate the higher QRPA 
effects). However, the interpretration of such calculations 
depends on the extra modelling, and the underlying ap- 
proximations, of collective excitations through the QRPA 
method; further studies are needed. 
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